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ABSTRACT 

Spectral line and continuum observations of the ionized and molecular gas in G20.08-0.14 N explore 
the dynamics of accretion over a range of spatial scales in this massive star-forming region. Very Large 
Array observations of NH3 at 4" angular resolution show a large-scale (0.5 pc) molecular accretion 
flow around and into a star cluster with three small, bright H II regions. Higher resolution (0.4") 
observations with the Submillimeter Array in hot core molecules (CH3CN, OCS, and SO2) and the 
VLA in NH3, show that the two brightest and smallest H II regions are themselves surrounded by 
smaller scale (0.05 pc) accretion flows. The axes of rotation of the large and small scale flows are 
aligned, and the timescale for the contraction of the cloud is short enough, 0.1 Myr, for the large-scale 
accretion flow to deliver significant mass to the smaller scales within the star formation timescale. 
The flow structure appears to be continuous and hierarchical from larger to smaller scales. 

Millimeter radio recombination line (RRL) observations at 0.4" angular resolution indicate rotation 
and outflow of the ionized gas within the brightest H II region (A). The broad recombination lines and 
a continuum spectral energy distribution (SED) that rises continuously from cm to mm wavelengths, 
are both characteristic of the class of H II regions known as "broad recombination line objects". 
The SED indicates a density gradient inside this H II region, and the RRLs suggest supersonic flows. 
These observations are consistent with photoevaporation of the inner part of the rotationally flattened 
molecular accretion flow. 

We also report the serendipitous detection of a new NH3 (3,3) maser. 
Subject headings: H II regions — ISM: individual (G20.08-0.14) — masers — stars: formation 



1. INTRODUCTION 

Massive star-forming regions (MSFRs) with O stars are 
usually identified by a group of hypercompact (HC) H II 
or ultracompact (UC) H II regions found together, deeply 
embedded in a dense molecular cloud (IGaray fc Lizanol 
[LQfll ICIra7cTr^l27)0l iHoare et al.ll2007r i. That several 
H II regions are typically found within each star-forming 
region indicates that massive stars form together in small 
clusters. Furthermore, the infrared luminosity and radio 
continuum brightness of the individual H II regions sug- 
gest that some of them may themselves contain more 
than one massive star. Thus, the spatial structure of 
massive star-forming regions is clustered and hierarchi- 
cal: the star-forming regions contain a number of sepa- 
rate HC and UC H II regions, each of which may in turn 
contain a few massive stars. 

Low angular resolution, single-dish, molecular line sur- 
veys of MSFRs show evidence for lar ge scale contraction 
of the embedding molec ular clouds ( Wu fc Evansl [20031 : 
iKlaassen fc W ilson 2008). Higher angular resolution ob- 
servations of some of these regions identify velocity gra- 
dients consistent with rotation and inflow. In addition 
to the accretion flows seen on the large-scale (~ 0.3 — 1 
pc) of the embedding molecular cloud (G10.6-0.4: Ho & 
Haschick 1986, Keto et al. 1987a, Keto 1990; G29.96- 
0.02: Olmi et al. 2003), accretion flows are also seen on 
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smaller (< 0.1 pc) scales around individual HC and UC 
H II regions (G10.6-0.4: Keto et al. 1988, Sollins et al. 
2005a; W3(OH): Keto et al. 1987b, Keto et al. 1995; 
W51e2: Zhang & Ho 1997, Young et al. 1998; G28.20- 
0.05: Sollins et al. 2005b; G24.78+0.08: Beltran et al. 
2004, 2005, 2006, Galvan-Madrid et al. 2008; G29.96- 
0.02: Beuther et al. 2007). 

It is unclear how the flows on different length scales are 
related. In the case of G10.6-0.4, the cluster-scale accre- 
tion flow can be traced down from the largest cloud scale 
to the small scale of the brightest H II region, but it is not 
known whether this holds for other objects. For exam- 
ple, in a survey of MSFRs, selected on the basis of IRAS 
colors and specifically excluding those with H II regions, 
multiple bipolar molecular outflows (implying the pres- 
ence of accretion flow s) are s een in random orientations 
(jBeuther et alj|2002alfb( . l200l . The different orientations 
of these smaller-scale flows suggest separate, individual 
centers of collapse. This comparison raises the question 
whether a large-scale coherent flow is required for the 
formation of the most massive stars, O stars (M* > 20 
Mq) capable of producing bright H II regions, whereas 
B stars require only smaller scale flows. 

It is also unclear what happens in an accretion flow 
when the inflowing molecular gas reaches the boundary of 
an embedded H II region. Previous observations suggest 
that the H II regions in an MSFR that are surrounded 
by accretion flows, may be best understood as deriv- 
ing from th e cont i nuous ionization of the accretion flow 
(|K^[2002l l200l [2001 . rather than as a dynamically 
separate expanding bubble of ionized gas within the flow. 
Part of the ionized gas may continue to the central star 
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or stars and part escapes off the rotationally flattened 
accretion flow as a photoevapor ative outflow perpendic- 
ular to the plane of rotation ( Hollenb ach et al.l 119941: 
Yorke fc Weldll996HLizano et al.lll996HJohnstone et all 
19981 iLugo et al.ll2004f ). The outflow is accelerated to 
supersonic speeds by the density gradient maintained 
by the stellar gravity (|Ketol [20071 ). Because the extent 
of an ionized outflow is generally larger than the re- 
gion of ionized inflow, in most cases the outflow should 
be detected more easily than the inflow. H II regions 
classified as "broad recombina t ion line objects" (B RLO) 
(jJaffe fe Martm-Pintadol 119991 : ISewilo et al l 12004ft show 
steep density gradients and supersonic flows (|Keto etldl 
2008), consistent with photoevaporation and accelera- 
tion. It is not known whether all BRLO are associ- 
ated with accretion. If the accretion surrounding an O 
star cluster is continuous from the largest to the smallest 
scales, this must be the case. 

There are only a handful of radio recombination line 
(RRL) observations that spatially resolve the ionized 
flow within an HC H II region. Velocity gradients con- 
sistent with outflow and rotation in the ionized gas 
have been pre viously reported for W3 (OH) (jKeto et alJ 
Il99fih. W51e2 (IKeto fc Klaassenl [2008h. and G28.20-0.05 
(|SewiIo et al.l l2008h~ Observations of the very massive 
and spatially large G10.6-0.4 H II region made at the 
VLA in the highest possible ang ular resolution are ab le 
to map the inflowing ionized gas (jKeto fc Woodl l2006). 

In order to study the accretion dynamics over a range 
of scales in a MSFR, from the cluster scale down to the 
scale of individual HC H II regions and within the ionized 
gas, we set up a program of radio frequency molecular 
line, recombination line, and continuum observations at 
two telescopes and with several different angular resolu- 
tions. For this study we chose the massive star formation 
region G20.08-0. 14 North (hereafter G20.08N), identified 
by three UC and HC H II regions detec ted in the cm 
continuum bv lWood fe Churchwelll (fl989l) . The total lu- 
minosity of the region is L ~ 6.6 x 10 5 Lq for a distance 
of 12.3 kpc 4 

Previous observations suggest accretion in the 
G20.08N cluster. Molecular-line obser vations show 
dense gas embedd ing the H II regions (jTurnerl 119791 : 
iPlume et alJ Il99*2i ). Molecular masers, generally asso- 
ciated with ongoing massive-star formation, have been 
detected in a number of studies (OH: Ho et al. 1983; 
H 2 0: Hofner fe Churchwell 1996 ; and CH 3 O H: Walsh 
et al. 1998). iKlaassen fc Wilson! (f2007l [2008h observed 
large-scale inward motions consistent with an overall con- 
traction of the embedding molecular cloud. Those au- 
thors also observed SiO line profiles suggestive of mas- 
sive molecular outflows, further evidence for accretion 
and star for mation. The recombination line spectra show 
broad lines (|Garav et all 119851 : ISewilo et al1l2004f) . pre- 

4 Both near and far kinematic d istances have been re ported for 
G20.08N. The near value given by IDownes et ail [Ti980h (d 4.1 
kpc) is t he most commonl y quo ted in the previous literat ure. In 
contrast, Fish et al. (2003) and Anderson & Bania (2009) report 
that this region is at the far kinematic distance (d 12.3 kpc). We 
will assume the far distance throughout the rest of the paper. For 
reference, a scale of 0.5" corresponds to PS 6000 AU (0.03 pc). The 
total luminosity of the region was esti mated to be L ~ 7.3 X 10 4 Lp , 
assuming the near kinematic distance (Wood & Churchwell 1989). 
Correcting for the location at the far distance, the luminosity is 
L ~ 6.6 x 10 5 L Q . 



sumably due to large, organized motions in the ionized 
gas. However, the previous observations do not have the 
angular resolution and the range of spatial scales needed 
to confirm the presence of accretion flows and study them 
in detail. 

In this paper we report on several observations of 
G20.08N and discuss our findings. We confirm active 
accretion within the cluster. Furthermore, we find that 
the parsec-scale accretion flow fragments into smaller 
flows around the individual HC H II regions, and that 
the gas probably flows from the largest scale down to 
the smallest scale. This continuous and hierarchical ac- 
cretion may be necessary to supply enough mass to the 
small-scale flows to form O-type stars, in contrast to low- 
and intermediate-mass star-forming regions with stars no 
more massive than ~ 20 Mq, where isolated accretion 
flows around individual protostars may be sufficient. 

2. OBSERVATIONS 

2.1. SMA 

We observed G20.08N on 2006 June 25 a nd July 6 with 
the Submillimeter Array 5 (|Ho et alJ l2004) in its very ex- 
tended (VEX) configuration. Two sidebands covered the 
frequency ranges of 220.3 - 222.3 GHz and 230.3 - 232.3 
GHz with a spectral resolution of « 0.5 km s - . The 
H30a recombination line (vq — 231.9009 GHz) was posi- 
tioned in the upper sideband 6 . The observations sampled 
baseline lengths from « 50 to < 400 kA, sensitive to a 
range of spatial scales from ks 0.5" to w 4.1". 

The visibilities of each observation were separately 
calibrated using the SMA's data calibration program, 
MIR. Table 1 lists relevant information on the calibra- 
tors. We used quasars for the absolute amplitude scale 
as well as the time-dependent phase corrections and 
frequency-dependent bandpass corrections. Inspection 
of the quasar fluxes and comparison with their historical 
flux densities in the SMA database 7 suggest that flux cal- 
ibration is accurate to better than 20%. The calibrated 
data were exported to MIRIAD for further processing 
and imaging. 

There were enough line-free channels in the 2-GHz 
passband to subtract the continuum in the (u, v) domain. 
The line-free continuum was self-calibrated in phase, and 
the gain solutions were applied to the spectral line data. 
We list all the identified lines in Table 2. Figure 1 shows 
the continuum-free spectra across the entire sidebands at 
the position of the 1.3-mm peak. 

To improve the sensitivity, the data were smoothed 
to a spectral resolution of 2 km s - . The rms noise 
in our natural- weighted maps, made from the combined 
observations of both days, is ~ 2 mJy beam -1 for the 
single-sideband continuum and ~ 30 mJy beam -1 per 
channel (2 km s -1 wide) for the line data. 

2.2. VLA 

Spectral line observations of NH3 (J,K) = (2,2) and 
(3,3) at two different angular resolutions were obtained 

5 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

6 In "chunk" 20 of the SMA correlator setup 

7 http:/ /smal. sma.hawaii.edu/callist/callist. html 
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with the Very Large Array 8 in the D and BnA con- 
figurations (projects AS749, AS771, and AS785). Par- 
tial results fr om these observations were presented in 
iSollinsI (|2005h . All the observations, except the VLA- 
D (3,3), were done with a bandwidth of 3.125 MHz 
(« 39 km s - ) divided in 64 spectral channels, each 
0.6 km s -1 wide. The VLA-D (3,3) observation was 
done with the same bandwidth divided in 128 spectral 
channels, each 0.3 km s -1 wide. The bandwidth covers 
the main hyperfine line and one line from the innermost 
satellite pair 9 . The NH 3 data are presented at a spectral 
resolution of 0.6 km s - . The noise per channel in the 
final images was in the range of 1.0 — 1.5 mJy beam -1 . 

In addition to the molecular line observations, we ob- 
served the ionized gas in the H66a recombination line in 
the VLA B configuration. The correlator was set up to 
cover a bandwidth of 12.5 MHz (« 166 km s -1 ) divided 
in 64 channels of 2.6 km s -1 each. The rms noise per 
channel in the final image was ~ 1 mJy beam -1 . 

All three VLA data sets were calibrated using stan- 
dard procedures in the AIPS software. Tables 1 and 2 
summarize the relevant observational parameters. The 
continuum was constructed in the (u, v) domain from 
line-free channels and was then self-calibrated. The gain 
solutions from self-calibration were applied to the line 
data. 

3. RESULTS AND DISCUSSION 
3.1. The Continuum Emission 
3.1.1. Morphology 

Figure [5] shows the 1.3-cm continuum (contours) ob- 
tained from the VLA-BnA observations overlaid with 
the 1.3-mm continuum from the SMA-VEX data (color 
scale). At 1.3 cm we resolve the G20.08N system into 
the th ree components reported by IWood fc Churchwell 
( 1989). H II region A is the brightest, westernmost peak. 
H II region B is the slightly broader peak w 0.7" to the 
SE of A. H II region C is the more extended UC H II fur- 
ther to the SE. Its brightest, eastern rim is detected at 
the ~ 10 mJy beam -1 level in our 1.3-mm observations. 

The continuum of H II region A is unresolved at 1.3 
cm; at 1.3 mm it shows a core- halo morphology. The 1.3- 
mm core is unresolved (Gaussian fits yield a deconvolved 
size at half power FWHM < 0.4"). The low-intensity 
halo has a diameter of » 1" (see Fig. H). The H30a 
emission is confined to the unresolved core (§ I3.4[) and 
the warm molecular gas (§ 13. 3p coincides with the ex- 
tended continuum halo. This indicates that the unre- 
solved H II region A is surrounded by a dust cocoon. H 
II region B is barely resolved in the 1.3-cm VLA BnA 
map (deconvolved FWHM « 0.6"). The peak position 
of H II region A is identical at 1.3 cm and 1.3 mm: 
a(J2000) = 18 h 28 m 10?30, c5(J2000) = -11° 28' 47"8, 
within the positional uncertainty of the reference quasars 
(in the range 0.01" - 0.1"). 

3.1.2. SED and the Nature of the Millimeter Emission 



8 The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 

9 The NH 3 mol ecule is symmetric top with inversion, see 
IHo fc Townesl (1983) for details. 



Previous observations of the recombination lines at 2 
and 6 cm (jSewilo et alj 12004) put the H II regions in 
G20.08N in the class known as "broad recombination 
line objects" (although those observations did not have 
sufficient angular resolution to separate H II regions A, 
B, and C). The large widths in cm- wavelength recombi- 
nation lines are due to pressure broadening at high gas 
densities (> 10 5 cm -3 ) as well as unresolved supersonic 
motions (|Keto et alj|200l . These H II regions also have 
continuum spectral energy distributions (SEDs) that in- 
crease with frequency through the mm wavelengths, ev- 
idence for a steep density gradient in the ionized gas 
(iFranco et al J 12000 : lKetoll200l lAvalos et alJl20M [20091 : 
iKeto et al.ll2008D . Our new VLA and SMA observations, 
at 1.3 cm and 1.3 mm respectively, extend the SED to 
millimeter wavelengths. We find that the flux density of 
H II region A continues to rise from cm to mm wave- 
lengths (Fig. [2]), and we analyze this characteristic in 
detail below. The flux densities of H II region B at 1.3 
mm (k, 93 mJy) and 1.3 cm (« 202 mJy) imply a spec- 
tral index (a, where S v oc v a ) of a ~ —0.3, roughly 
consistent with the expected index of —0.1 of optically 
thin gas. H II region C is more extended and most of the 
1.3-mm flux is resolved out. 

At wavelengths shorter than 1 mm, thermal dust emis- 
sion contributes significantly to the continuum. We esti- 
mate the relative contributions of dust and free- free emis- 
sion at 1.3 mm from our recombination line observations 
(described in § I3.4j) and the theoretically expected line- 
to-continuum ratio. In the optically thin limit (a good 
approximation at 1.3 mm) the free-free linc-to-continuum 
flux rat io Sl/Sq is given by the ra tio of the opacities 
h l /k c ([Gordon fc Sorochenkoll2002ft : 



kl 



cm 



x exp 



Kc 



= 9.77 x 10 



_3 T^e^z 



(27rm,/fc)VW™' TV* 
'§*)(! -.-./«■), (la) 

T 3/2 



2 T 3/2 



17.72 + In- 



(16) 



where all the units are in cgs, the physical constants 
have their usual meanings, n\ = 30 for H 30a, f ni ,n 2 ~ 
0.1907ni(l + 1.5/ni) for a lines, n e = ni for hydrogen, 
T w 8000 - 10 000 K, and cf> v is the normalized line pro- 
file. The main source of uncertainty in equations [1] and 
[T] is the temperature of the ionized gas. Assuming that 
the line profile is Gaussian and correcting for 8 % he- 
lium in the gas, the expected ratio at the line center is 
S L o/S c ~ 3.0 for T = 8000 K, or S L /S C « 2.3 for 
T = 10 000 K. The observed ratio is 'S Lfi /S c ~ 1.8. 
Therefore, assuming that the RRLs are in LTE, the free- 
free contribution to the 1.3-mm flux of H II region A is 
w 60% — 80% for the assumed temperature range. 

Figure [3] shows a model SED for H II region A in which 
70 % (355 mJy) of the 1.3-mm flux is produced by free- 
free from an H II region with a density gradient and 30 
% (142 mJy) by warm dust. Assuming radiative equi- 
librium, we set the dust temperature Td to 230 K, the 
average temperature of the dense gas at the same scales 
(i j3.3.2p . The modeling procedure is described in iKetol 
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(2003) and IKeto et all ([2008h . Table 3 summarizes the 
model. The total gas mass inferred from the dust emis- 
sion is too large (M ^35 — 95 Af Q ) for the H II region 
alone, so most of the dust must be in the cocoon around 
the H II region. The calculated mass range takes into 
account uncertainties in the dust emissivity, but not in 
the temperature. 

The density gradient derived for the ionized gas in H 
II region A is n e oc r~ 7 , with 7 — 1.3. Equilibrium 
between recombination and ionization in this model H II 
region requires an ionizing fl ux equivalent to a n 07.5 star 
(using the computations of iVacca et al.lll996f ). although 
this ionizing flux could be made up of several stars of 
slightly later spectral type. The model for the SED as- 
sumes spherical geometry and a static H II region with no 
inflow of neutral gas into the H II region. In contrast, H 
II region A is embedded in a rotationally flattened accre- 
tion flow (see § 13. 3p . so the determination of the stellar 
spectral type is only approximate. 

3.2. The Large- Scale Molecular Cloud 

The large-scale molecular cloud is detected by the NH3 
VLA-D observations. From these data we find the pres- 
ence of a parsec-scale accretion flow surrounding the clus- 
ter of H II regions. In § 13.2.11 we first estimate the 
systemic velocity V sys of the cloud with respect to the 
local standard of rest (LSR), and then analyze the line 
velocities to determine rotation and infall following the 
procedure used in mea suring accretion flow yelocities in 



previous papers (e.g.. iHo fc Haschickl 
ll987aHZhang fc Hdfl997l lYoung et alj 
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IKeto et al.1 
. In § I3~2T21 

we derive the physical properties of the cloud: tempera- 
ture, mass, ammonia abundance, and density. 

3.2.1. Dynamics 

Figure [4] shows the channel maps of the (J, K) = (3, 3) 
main hyperfine line. The most notable feature is a veloc- 
ity gradient along the major axis of the cloud, consistent 
with rotation. The redshifted emission is toward the NE, 
while the blueshifted gas is toward the SW. Also, there 
is strong absorption against the free-free background of 
the H II regions at the center of the cloud, which are un- 
resolved at ~ 4" resolution. The separation between the 
emission peaks on either side of the absorption is ks 10", 
or 0.6 pc. The symmetry of the channel maps suggests 
a systemic velocity of about 42 km s _1 . A Gaussian fit 
to the emission spectrum integrated in a box covering 
the entire cloud, and clipping out the redshifted absorp- 
tion, gives a centroid velocity of 42.3 ± 0.4 km s" 1 . We 
adopt a systemic velocity for the parsec-scale cloud of 
V sys — 42 km s _1 , in agreement with that reported by 
iPlume et al.l |l992) (41.9 km s _1 ) based on observations 
of CS J = 7 - 6. 

Figure [5] shows position- velocity (PV) diagrams across 
cuts at PA = 45° (SW-NE, major axis) and PA = 135° 
(NW-SE, minor axis). The rotation is seen in the SW- 
NE cut as a shift in the velocity of the emission contours 
from one side of the absorption to the other. The velocity 
offset with respect to V sys seen in the emission contours 
increases inward, suggesting that the gas rotates faster 
with decreasing radius. 

Under the assumption that the velocity gradient seen 
in emission along the major axis is dominated by rota- 



tion, the redshifted absorption in the PV diagrams (Fig. 
[5]) is also evidence for inward flow toward the central H 
II regions, as there is an excess of redshifted absorption. 
This is more clearly seen in Fig. [6j which shows the spec- 
tra toward the absorption center. The NH3 (2,2) main 
hyperfine absorption peak is redshifted by 2.3 kms -1 
with respect to the systemic velocity. An infall veloc- 
ity Vi n f « 2 km s _1 is also seen in the NH3 (3,3) line, 
although the spectrum is contaminated by an NH3 (3,3) 
maser (see description in Fig. [6]). The maser is confirmed 
by our high angular resolution NH3 observations (§ 13. 6| . 

If the rotation were seen edge-on, there would not be a 
velocity gradient across the H II region along the minor 
axis (NW-SE). At an oblique viewing angle a velocity 
gradient along the minor axis is created by the inflow. 
The velocity gradient along both the minor and major 
axes (Fig. 5) implies that the rotationally flattened flow 
is tipped with respect to the line of sight. 

3.2.2. Cloud Parameters 

The optical depth of the gas can be determined from 
the brightness ratios of the hyperfine lines. From their 
optical depth ratio the rotational temperature between 
the (2, 2) and (3, 3) transitions c an be determined (see 
IKeto et al.lfl987at iMangum et 11111991 for details of the 
procedure). This temperature can be considered as a 
first-order approximati on to the kine t ic tem perature (Tfe) 
of the molecular gas. iDanbv et al.l (|1988f ) suggest that 
an improved estimate of the kinetic temperature (T^) is 
obtained by accounting for the populations in the up- 
per states of the K-ladders rather than assuming that all 
states above the lowest are negligibly populated. Includ- 
ing this correction, we obtain an average kinetic temper- 
ature Xfc ~ 22 K for the parsec-scale cloud. The high- 
resolution data discussed in § 13.31 show that higher gas 
temperatures are found closer to the H II regions. 

The mean column density of NH3 is Nnh 3 ~ 2.3 x 10 16 
cm~ 2 . We can determine the abundance of ammonia 
X^NHs) by comparing the NH3 and H2 column densi- 
ties. The mass of the cloud can be estimated from the ob- 
served velocity dispersion (oy ~ 3.5 kms -1 ), the radius 
of the cloud (R ~ 100 000 AU), and the virial theorem: 
M mr = (5/3G)i?c4 ~ 2300 M Q . Similarly, from the ob- 
served velocity gradient, the rotation velocity is V ro t ~ 4 
km s -1 at a radius R ~ 100 000 AU. Equating the cen- 
tripetal and gravitational forces, the gas mass inside R 
is M = {RV? ot )/G ~ 1800 M Q . Therefore, we estimate 
the average column density of H2 to be Njj 2 ~ 5 x 10 22 
cm~ 2 and X(NH^) ~ 5 x 10~ 7 , withi n a factor o f three 
of the value estimated for G10.6-0.4 (|Ketolll990D . The 
mean H2 density in the large-scale cloud is uh 2 ~ 10 4 
cm~ 3 . 

3.3. Molecular Gas in The Inner 0.1 pc 

Our SMA Very Extended (VEX) and VLA-BnA obser- 
vations provide a rich view of the molecular environment 
close to the H II regions. The lines detected with the 
SMA that have peak intensities > 300 mJy beam -1 and 
are not blended with any other line are listed in the top 
part of Table 2. We follow a similar outline as for the 
large-scale cloud. 

3.3.1. Dynamics 
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Figure [7] shows the channel maps of OCS J = 19 — 18, 
S0 2 J(K a , K b ) = 11(1, 11) - 10(0, 10), CH 3 CN J{K) = 
12(3) - 11(3), and CH 3 CN J(K) = 12(4) - 11(4). These 
maps show that the line emission from these hot-core 
molecules is considerably brighter around H II region A. 
All the molecules show a velocity gradient across this 
source, from the southwest to northeast, in a similar ori- 
entation as the larger, cluster-scale flow. We do not de- 
tect molecular emission around H II region B. At this 
high angular resolution we are not sensitive to bright- 
ness temperatures of less than ~ 10 K (3ct) for the mm 
lines. 

Figure [5] shows the velocity-integrated intensity (mo- 
ment 0) and the intensity-weighted mean velocity (mo- 
ment 1) maps of the four lines previously shown in Fig. 
[3 The integrated emission of CH 3 CN 12(4) - 11(4) is 
brightest in front of the H II region, while for the other 
molecules the brightness peak is slightly offset from the 
continuum. These differences reflect the relative bright- 
ness of each molecule with respect to the continuum emis- 
sion of H II region A. The velocity gradient across H 
II region A is also seen. Figure [9] shows the velocity- 
integrated (moment 0) and velocity dispersion (moment 
2) maps for the same lines as Fig. [8l The line widths 
increase toward the continuum peak, indicating that un- 
resolved motions increase closer to the H II region. 

Figure 10 shows the position-velocity (PV) diagrams 
for the lines of Figure 7 in cuts at position angles PA 
= 45° and PA = 135° across the continuum peak of HII 
region A. The cuts at 45° show the velocity gradient also 
seen in the channel maps of OCS and both of the CH3CN 
lines. The PV diagram of CH 3 CN 12(4)-11(4) has a fea- 
ture suggestive of a velocity gradient in the perpendicular 
direction PA = 135°, with an excess of redshifted emis- 
sion toward the NW. Consistent with our interpretation 
of the NH 3 VLA-D data, this suggests inward motion in a 
rotationally-flattened flow that is seen not quite edge-on. 
The same redshifted feature is also marginally detected 
in the lower excitation CH3CN transition as well as the 
OCS line. However, the infall signature in emission is 
only tentative, and a clearer indication of infall at small 
scales comes from the redshifted NH3 absorption in the 
VLA-BnA data (see below). 

In general, observational experience suggests that 
CH3CN, along with NH3, is a reliable tracer of 
high-density m ol ecula r mass and ac cretio n flows 
(ICesaroni et all fl997l: IZhang et all fl998l 12001 
iPatel et all 120051 ). CH 3 CN has recently been de- 
tected in the outflow of the nearby low-mass star 
LI 157 (|Arce et al.l [2008h . but at a very low brightness 
(< 0.03 K). The distribution of the OCS molecule in 
our observations is very similar to the CH 3 CN, but 
the SO2 velocities do not show the same pattern, and 
are more difficult to interpret. The SO2 may be more 
easily affected by the excitation conditions, and part 
of the observed emission could arise from the shocked 
boundaries of outflows. From Gaussian fits to the 
CH3CN K — 2,3, and 4 emission lines at the position 
of the continuum peak, the systemic velocity at small 
scales is estimated to be V sys = 41.8 ± 0.3 km s~ . 
Figure [TT] shows the CH3CN spectra and their Gaussian 
fits. 

At subarcsecond angular resolution, our NH 3 obser- 



vations are sensitive to emission of brightness tempera- 
ture above ~ 200 K. Therefore the thermal NH 3 is de- 
tected only in absorption against the bright continuum. 
As in § 13.2.21 a comparison between the NH3 absorp- 
tion line velocity and V sys shows an inward velocity of 
w 2 km s _1 in front of H II region A (Figure [12]) . The 
NH3 (2,2) absorption line in front of H II region B is red- 
shifted by ~ 2.0 km s - with respect to V sys , implying 
inward motion and accretion toward H II region B as well 
(Figure IT2")). The (3,3) absorption in front of H II region 
B is mixed with NH3 maser emission (§ 13. 6p . and the 
determination of the inward velocity is uncertain. More 
sensitive observations are needed to constrain the prop- 
erties of the molecular gas around H II region B. 

3.3.2. Core Parameters 

We estimate the dynamical mass M (gas plus stars) 
within the smaller accretion flow in the same way as with 
the large-scale flow (see § I3.2.2|) . At a radius R ~ 5000 
AU the rotation velocity is V ro t ^3 — 4 km s - . There- 
fore, M ~ 50 — 90 Mq . This is consistent with the lower 
limit to the stellar mass, M* « 35 Mq, required for 
ionization equilibrium (§ I3.1.2p . The estimate is also 
consistent with the gas mass derived from the mm con- 
tinuum once the free-free contribution has been properly 
subtracted, M gas ~ 35 - 95 Mq ft J3~l~2]) . The mean H 2 
density is nn 2 ~ 10 6 cm~ 3 . 

We derive the temperature in the dense gas surround- 
ing H II region A from the rotational energy diagram 
of t he lines of the CH3CN J = 12 — 11 K-lad ders (see 
e.g., lLoren fc Mundvl 119841: IZhang et all I1998D . Figure 
[T3l shows this diagram for two cases: one considering 
all the K = 0, ...,7 lines, and the other including only 
the K — 4,..., 7 lines, which have lower optical depths 
than the low-number K lines (also, the K = 0, 1 lines are 

blended, separated by only 5.8 km s _1 ). The rotational 
temperature obtained for the former case is T rot ~ 403 
K, while for the latter it is T rot ~ 230 K. The differ- 
ence between the two values given above appears to be 
caused by optical depth effects (the rotational diagram 
analysis assumes that the emission is optically thin) . Al- 
though our sensitivity level does not permit us to detect 
the lines of the isotopologue CH 13 CN and measure the 
optical depth of the CH3CN emission, the upper limits 
are not restrictive (tch 3 cn < 20). High optical depths 
are also suggested by the flat slope of the K ladders (Fig. 
1). F itting the K ladders taking into account the opac- 
ities (Qiu & Zhang 2009) yields a kinetic temperature 
close to the lower estimation, T% ~ 225 K, and optical 
depths are >> 1 (Keping Qiu 2009, personal communi- 
cation). The rotational temperatures are higher if we use 
only the brightest pixels instead of averaging over all the 
emission. This suggests a temperature gradient toward 
the center of the H II region. 

We derive a CH 3 CN column density of Nch 3 cn ~ 
7.5 x 10 15 cm -2 assuming Tk ~ 230 K. Comparison of the 
CH3CN column density with the dynamical mass implies 
an abundance A(CH 3 CN) - 5 x 10~ 9 - 2 x 10~ 8 for the 
range of masses quoted above. Abundance estimates in 
other MSFRs cover a range of values: ~ 10~ 10 inside the 
Orion hot core and ~ 10~ u outside (jLoren fc Mundvl 
119841) : 10~ 8 inside the Orion hot core and 10" 9 - 10" 10 
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in the Orion ridge (|Wilner et all 11994); - 3 x 10~ 8 in 



Sgr B2(N) (INummelin et al.l [2000h ; - 1 - 5 x 10" 7 in 
W51el/e2 (jRemiian et al.H2004h . 

The optical depth of NH3 in absorption toward H II 
region A is T2.2 ~ 4.2, t^^ ~ 0.9. The mean col- 
umn density is Nnh 3 



3.6 x 10 16 cm" 2 . 



If the am- 



monia abundance at these small scales is in the range 
X(NH 3 ) = 2 x 1(T 6 - 1 x 1CT 7 , then the molecular hy- 
drogen column density is N H2 ~ 2x 10 22 — 4 x 10 23 cnT 



This implies a molecular gas mass of M gas 



40 

Mq. The kinetic temperature obtained from the NH 3 
is Tfc ~ 50 K, considerably cooler than that obtained 
for CH3CN and implying that most of the NH3 column 
density is further away from the H II region than the 
CH 3 CN. 

The distribution of the warm molecular gas around H 
II region A appears to be rotationally flattened (Figs. El 
EJ, but the observed size is slightly larger than the 
radius at which the accretion flow is expected to become 
centrifugally supported (i.e., a "disk"), Rd = GM/V^ ot ~ 
2000 - 3000 AU for a star of mass M* = 35 M Q and 
V ro t = 3 — 4 km s . Actually, the disk scale matches 
the size of the H II region (Table 3, i J3~4|) . 

The mass-inflow rate toward H II region A can be es- 
timated from the high-resolution NH3 absorption. From 
an inflow velocity of Vi n j ~ 2 km s _1 and spherical ge- 
ometry, the mass-inflow rate is M ~ 1 x 10~ 3 — 2 x 
10 -2 M Q yr^ 1 (for the 2-40 M Q of molecular gas de- 
tected in ammonia absorption) . This estimate may be an 
upper limit because the hot molecular core is flattened 
rather than spherical. 

3.4. The Ionized Gas: Radio Recombination Lines 



In § 13.1.21 we inferred a density gradient inside the 
HC H II region A. In this section, we derive the inter- 
nal dynamics of this H II region based on multifrequency 
RRLs. The mm/sub-mm lines are especially important 
because they are much less affected by pressure broaden- 
ing and preferentially trace denser gas. While subarcsec- 
ond resolution studies at wavelength s longer than 7 mm 
have been available for many years (|De Pree et al.llT997l . 
12004 iRodriguez et al.ll200"9l ). similar studies at shorter 
wavelengths have had limited angular reso l ution (e.g., 
Uaffe fe Martin-Pintado! 11999ft . iKeto et alj (|2008l ) pre- 
sented the first results of high-frequency, high-resolution 
(~ 1"), multifrequency RRL observations in a sample 
of 5 MSFRs with similar characteristics to G20.08N A. 
They were able to separate the contributions of pressure 
broadening and large-scale motions to the line width, 
even when the H II regions were unresolved. We follow 
their procedure to analyze our RRL data. 

We observed the H30a (v = 231.90 GHz) and H66a 
(vq = 22.36 GHz) lines at subarcsecond angular resolu- 
tion (see Table 2). Because both the line-to-continuum 
ratio and the continuum intensity are lower at 22 GHz 
than at 231 GHz, the H66a line is much weaker than 
the H30a line. This is somewhat alleviated by the bet- 
ter sensitivity of the VLA, but the signal-to-noise ratio 
(S/N)in the high-frequency line is still better. Figure fT4l 
shows the moment and moment 1 maps of the H30a 
line. Although the emission is unresolved at half power, 
there is a slight indication of a velocity gradient in the 
ionized gas that agrees (not perfectly) with the rotation 



seen in CH 3 CN and OCS. Figure □>] shows the H30a 
and H66a spectra toward H II region A and their Gaus- 
sian fits. The H66a line shows evidence of a blueshifted 
wing, suggesting either inflow or outflow in addition to 
rotation. Within the uncertainties of the fits, both the 
H66a and H30a lines have the same line width (Table 
4). 

Assuming that the dynamical broadening Avd (caused 
by turbulence and ordered motions) and the thermal 
broadening Avt are Gaussian, and that the pressure 
broadening Avl is Lorentzian, the RRL has a Vo igt pro- 
file with line width ([Gordon fe Sorochenkoll2002ft : 



Av v {u) w 0.534Au L (i/) + (Av 2 D 



A4 + 0.217Av?») 1/2 , 



(2) 



where all the widths are FWHM. 

For the H30a line at 231.9 GHz the pressure broad- 
ening is less than the thermal broadening at e lectron 
densities below 2 x 10 8 cm" 3 (jKeto et al.l 12008ft . Our 
SED modeling indicates lower densities over most of the 
HII region. Therefore, the observed line width can be 
attributed to thermal plus dynamical broadening. The 
electron temperature T e in UC H II regions is typically 
T e = 8000 - 10 000 K, with a small gradient as a func- 
tion of Galactocentric radius (A fflerbach et al.lll996l ) . We 
adopt T e = 9000 K (T B < 100 K for H30a because of the 
low optical depth and < 1 filling factor), which trans- 
lates into a thermal FWHM of Av T = 20.9 km s" 1 . 
Therefore, from Eq. [TJ we obtain a dynamical width of 
Avd — 24.8 km s _1 . From the velocity gradient (Fig. 
[LT| it is seen that ~ 6 km s _1 of Avd can be in the 
form of rotation. The rest could be caused by inflowing 
or outflowing ionized gas, as suggested by the blueshifted 
(^2 — 3 km s _1 ) mean velocities of the RRLs, and by 
the blue wing in the H66a spectrum (Fig. I15[) . 

Most of the ionized gas that we see should be out- 
flow. Inflow inside the H II region is expected within the 
radius where the escape velocity from the star exceeds 
the sound speed of the ionized gas. Thi s is approx i- 
mately the Bondi-Parker transonic radius (jKeto! [2007). 
R b = GAhjlcl = 5.5 AU M*/M , or about 190 AU 
for H II region A, assuming a sound speed c s = 9 kms -1 
and stellar mass M+ = 35 Mq . H II region A extends out 
to ~ 2500 AU, so most of the gas is not gravitationally 
bound to the star and flows outward. In this model, the 
outflow is continuously supplied by photoevaporation off 
the rotationally flattened accretion flow (T. Peters et al. 
2009, in preparation). The somewhat misaligned velocity 
gradient in the ionized gas (Fig. 14) derives from a com- 
bination of the rotation and outflow blended together in 
the observing beam. 

3.5. Outflow Tracers 

Low angular resolution (HPBW = 14"), single- 
dish (JCMT) observations of SiO (8-7) show evi- 
dence for large-scale mo lecular outflows in G20.08N 
(|Klaassen fc Wilson! [2008D . Although the standard out- 
flow tracers 12 CO J = 2 - 1 and 13 CO J = 2 - 1 are 
in our passband, we do not detect any CO in emission. 
The (it, v) coverage of our SMA-VEX observations is in- 
capable of imaging structures larger than ~ 4"; therefore, 
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the CO emission from the molecular cloud and molecu- 
lar outflow must be on larger scales. 12 CO and 13 CO 
are seen in our data in absorption at the position of H 
II region A, at several different velocities in the range 
Vlsr = 42 — 84 km s _1 . Some CO absorption features 
are at the same v elocities as the HI absorption features of 
iFish et all (|200l and are therefore due to foreground gas 
that is not related to G20.08N but rather to intervening 
Galactic spiral arms. 

3.6. A New NH 3 (3,3) Maser 

A handful of NH3 masers have been reported in the 
literature, always associated with massive star forma- 
tion (e.g., IWalsh et ani2007f) . Many of the known NH 3 
masers are from non-metastable (J > K) transitions. 
The first clear detection of a metastable (J = K) NH3 
(3,3) maser was reported by Mangu m fe Woottenl (|1994T ) 
toward DR 21 (OH). Most of the detections point to- 
ward a shock excitation origin for the population inver- 
sion, inasmuch as the maser spots are invariably associ- 
ated with outflow indicators such as bipolar CO and/or 
SiO str uctures, class I methanol masers, and/or water 
masers dMangum fc Woottenl 1199*41 : iKraemer fe Jacksonl 
[19951 IZhang fe Hoi 1 1995ft ~ 

We report the serendipitous detection of a new NH3 
(3,3) maser toward G20.08N. The maser spot is rela- 
tively weak, and is spatially centered at a(J2000) = 
18 h 28 m 10?346, (5(J2000) = -11° 28' 47"93, close in 
projection to the center of H II region B. The maser spot 
is spatially unresolved even in uniform weighting maps 
of the VLA-BnA data (HPBW = 0.48" x 0.26", PA = 
—0.8°). If the deconvolved source size is limited to half 
the beam size, then the peak brightness temperature of 
the spot is constrained to Tb > 7 x 10 3 K. The high 
intensity, together with the absence of similar emission 
in our (2,2) maps at high angular resolution, confirm the 
maser nature of the (3,3) emission. The spectral feature 
is also very narrow (Fig. 1 16[) . typical of maser emission, 
although it shows evidence of line wings. From a Gaus- 
sian fit to the line profile, the velocity of the maser is 
V m aser ~ 44.7 km s -1 (accurate only to 20 % - 40 %). 
The FWHM is ~ 0.7 km s _1 , after deconvolving the 
channel width of 0.6 km s _1 . Owing to its position, it is 
probable that the maser is excited by H II region B. We 
do not have sufficient data, however, to assert that it is 
excited in a shock. 

4. HIERARCHICAL ACCRETION IN G20.08N 

4.1. The Observations 

One of the central questions in star formation is 
whether star formation is "bimodal", i.e., whether 
high-mass and l ow- mass stars form in a different way 
(|Shu et al.l 119871 ). Analysis of recent observations sug- 
gests that accretion flows around protostars of all masses 
can be explained within the standard model of star for- 
mation as a combination of a thin disk inside a rota- 
tionally flattened envelope (Kumar & Grave 2007; Moli- 
nari et al. 2008; E. Keto & Q. Zhang 2009, in prepa- 
ration). Of course, once a protostar gains the mass and 
temperature of an O star, the formation of an H II re- 
gion within the accretion flow introduces new phenomena 
(|Keto|[2002l: [Keto fc Woo1[2^01IK^[2007T) . Yet not all 



molecular clouds form ma s sive s tars. In a comprehen- 
sive survey. [Solomon et alJ (1987) found that clouds that 
form low-mass stars are uniformly distributed through- 
out the Galactic disk, but those that form the clusters 
of the most massive stars, O stars capable of producing 
significant H II regions, are associated with the galactic 
spiral arms. If the formation process is very similar (with 
the addition of significant ionization for M* > 20M©), 
then the difference may be found in the conditions in 
the molecular clouds. The observations of G20.08 N re- 
ported here and of G10.6 -0.4 reported previously (jKetol 
fl990tlKeto fc Wood! [20061 ) suggest that one difference is 
that the molecular clouds surrounding young clusters of 
O stars are in a state of overall collapse whereas in star- 
forming regions without H II regions we see only localized 
collapse. 

In G20.08N the H II regions A, B, and C are sur- 
rounded by a common molecular cloud of radius ~ 0.5 
pc and mass M ~ 2000 M Q , which is rotating with a 
velocity of ~ 4 km s~ and contracting with an inward 
velocity of ~ 2 km s~ . This rotation and contraction 
constitute a large-scale accretion flow. The velocity of 
the inward flow is about equal to the rotational velocity 
implying that the gas is approximately in free-fall and 
not constrained by centrifugal force. Within this larger 
flow are at least two smaller accretion flows around H II 
regions A and B. The molecular core around H II region 
A is bright enough to be detected in hot-core molecules in 
emission at subarcsecond angular resolution. The radius 
and mass of this core are 0.05 pc and ~ 20 M Q . Accre- 
tion in the core is indicated by rotation at a velocity of 3 
or 4 km s _1 and contraction of ~ 2 km s~ . We do not 
detect molecular emission around H II region B, but the 
accretion inflow is inferred from NH3 absorption that is 
redshifted by ~ 2 km s _1 with respect to the systemic 
velocity. 

In contrast, star- forming regions that contain only 
low mass stars or even stars as massive as type B 
do not appear to have this global collapse of the en- 
tire parent cloud. The observations of these regions 
suggest only localized accretion flows within individ- 
ual clumps. For e xample, star-forming regions such as 
IRAS 19410+ 2336 ([Sridharan et al.ll2002t Beuther et al l 
20021 1 [2001 IRA S 05358+3543 d Beuther et alJT2 002a: 
Leurini et alTl2007l) . and AFGL 5142 ([Zhang et alj|200a . 
2007ft . have bolometric luminosities of at most a few 
times 10 4 Lq, consistent with type B protostars. There 
are a number of cores of size similar to those in G20.08N, 
but there is no reported evidence of a larger parsec-scale, 
accretion flow. 

In an analysis of a recent numerical simulation, 
iVazquez-Semadeni et alJ (|2009ft also find that the for- 
mation of massive stars or clusters is associated with 
large-scale collapse involving thousands of M© and ac- 
cretion rates of 10 -3 M© yr _1 . In contrast, low- and 
intermediate-mass stars or clusters in their simulation 
are associated with isolated accretion flows that are a 
factor of 10 smaller in size, mass, and accretion rate. 

If global collapse of the host molecular cloud is nec- 
essary for the formation of O stars in clusters, then the 
association of O stars with galactic spiral arms may im- 
ply that compression of giant molecular clouds as they 
pass through galactic spiral arms may be the primary 
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mechanism for i nitiating global collaps e (Roberts 119691 : 
iShu et all [1973 iMartmez-Garcia etall 120091). Th e low- 
mass star-forming regions found bvlSolo mon etaD ((1987? ) 
to be spread throughout the Galaxy may not need such 
large-scale compression. 

4.2. Resupply of the Star-Forming Cores 

The orientations of the large cloud-scale accretion flow 
and the core-scale flow around H II region A are similar 
and the flows could be continuous (we do not know the 
orientation of the flow around H II region B because we 
cannot detect the surrounding molecular emission). The 
molecular core around H II region A contains only a few 
tens of Mq, similar to the mass of a single O star. If 
the core is to form one or more O stars at less than 100 
% efficiency, its mass must be resupplied by the larger 
scale accretion flow. Resupply is also suggested by the 
short dynamical or crossing timescale, given by the ratio 
of the size to the infall velocity, ~ 10 4 yrs. If the core is 
to last more than this, it must be resupplied by the larger 
scale flow, which has a dynamical timescale of ~ 10 5 yrs. 
If the core is not resupplied then the growing protostar 
may simply run out of gas before reaching the mass of 
an O star. 

A similar process of resupply is suggested in recent 
theoretical work. Analytic arguments show that as an 
unstable cloud fragments, there should be a continu- 
ous cascade of mass from larger to s maller fragments 
as well as a cascade of kinetic energy (jField et all 12008: 
iNewman fc Wasserman|[l990[ ). In re cent numerical sim - 
ulations of high mass star formation. iWang et al.l (|2009l ) 
find that most of the mass is s upplied from outside a 0.1 
pc cor e around the protostar. IVazquez-Semadeni et ahl 
(2009) and T. Peters et al. (2009, in preparation) also 
find in their simulations that as a massive core is con- 
sumed from the inside by an accreting protostar, the core 
continues to accrete mass from the outside. 

4.3. Accretion Rate 

The accretion rate within the small-scale flow around 
H II region A is M ~ 1 X 10~ 3 to 2 x 10~ 2 M Q . Similar 
accretion rates are reported for flows around other H I I 
regions such as G10.6-0.4, 8x 10~ 4 M© yr" 1 (jKetoll 19901), 
and W51e2, 3 x 10~ 3 to 1 x 10~ 2 M© yr" 1 (|Young et all 
1998). In contrast, the accretion rates estimated for the 
cores in these MSFRs without bright H II regions are 
generally lower by one or two order s of magnitude, ~ 
10" 4 Mq yr" 1 (|Beuther et aJJl2002allbh . 

A high accretion rate is necessary to form an O star. 
Because accreting massive protostars begin core hydro- 
gen burning well before reaching the mass of an O star, 
they evolve essentially as main-sequence stars of equiva- 
lent mass. Numerical simulations of stellar evolution that 
include accretion predict that unless the accretion rate 
is high enough, a growing massive protostar will evolve 
off the main-seque nce and explode bef ore it reaches the 
mass of an O star (|Keto fc Woodl[20M ). For a protostar 
to gain the mass of an O star, the rate at which ac- 
cretion supplies fresh hydrogen to the growing protostar 
must be greater than the rate at which the star burns 
the hydrogen. At the upper end of the mass spectrum, 
M* > 40 M Q , this rate is > 10~ 3 M© yr" 1 . The ac- 
cretion flows in G20.08N are capable of supplying gas at 



the rate necessary to form massive O stars. In MSFRs 
without O stars, the accretion rates may be too low for 
the stars to achieve the mass of an O star within their 
hydrogen-burning lifetimes. 

4.4. Transfer of Angular Momentum 

If the flow in G20.08N is continuous, then the observa- 
tions show that the flow spins up as it contracts. Ignoring 
projection effects, the magnitude of the specific angular 
momentum on the large scale is L/m ~ 0.6 km s _1 pc 
atrsj 0.2 pc, from the VLA-D NH3 data, while from the 
VEX SMA data it is L/m ~ 0.1 km s" 1 pc at r « 0.02 
pc around H II region A. From these estimates ~ 85 % of 
the specific angular momentum in the large-scale flow is 
lost. In previous observations of G10.6-0.4 (Keto 1990), 
we found that 97% of the angular momentum in that flow 
is lost between 1.5 and 0.02 pc. Evidently, angular mo- 
mentum is efficiently transferred outward and does not 
prevent collapse of the cloud. 

Magnetic fields may be important in this process (e.g., 
iGirart et al.ll2009l) . Although there are no observations 
of the magnetic field in G20.08N, the field direction has 
been mapped in another MSFR with HC H II regions, 
W51el/e2. Here, observations show that the magnetic 
field is uniform on the larger scale, 0.5 pc (|Lai et al.l 
120011) . while on the smaller scale, 0.03 pc, of the accretion 
flow onto W51e2, the field is pinched into a n hourglass 
shape with the accretion flow at the waist (|Tang et al.l 
2009). Thus ordered field lines extend from the H II 
region-scale accretion flow to the large-scale molecular 
cloud, and if the field has enough strength, angular mo- 
mentum could be transferred outward by the field. How- 
ever, dust polarization observations do not give a di- 
rect estimate of the field strength, and it is also pos- 
sible that the field is essentially passive and just carried 
along by the flow. The observation that the clouds in 
both W51el/e2 and G20.08N are close to free-fall col- 
lapse implies that the magnetic field is not strong enough 
to support the clouds. In other words, these clouds are 
magnetically supercritical. 

Numerical simulations of star formation that do not 
include magnetic forces show that angular momentum 
can b e transferred by hydrodynamics alone. lAbel et alJ 
(2002) find that at any radius, there is both low and 
high angular momentum gas, and that pressure forces 
or shock waves can redistribute the angular momentum 
between fluid elements. Lower angular momentum gas 
sinks inward and displaces higher angular momentum gas 
outward, resulting in a net outward flow of angular mo- 
mentum. 

Whether the specific angular momentum is transferred 
by hydrodynamics or magnetic forces, observations and 
simulations show that angular momentum is not con- 
served as a function of radius, and does not prevent the 
gas from flowing continuously from large to small scales 
in a rotating flow. 

5. CONCLUSIONS 

We report radio and mm observations of the molecular 
and ionized gas toward the O-star cluster G20.08N, made 
with an angular resolution from ~ 0.1 pc to ~ 0.01 pc. 
Our main findings can be summarized as follows: 

1. We find a large-scale (~ 0.5 pc) accretion flow 
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around and into a star cluster with several O-type 
stars, identified by one UC and two HC H II re- 
gions. This flow is rotating and infalling towards 
its center. The two HC H II regions are surrounded 
by smaller accretion flows (~ 0.05 pc), each of them 
with the signature of infall too. The brightest (to- 
ward H II region A) is detected in mm emission 
lines, and rotates in concordance with the large- 
scale flow. 

2. The similar orientations of the flows at small and 
large scales, as well as their dynamical timescales 
(<~ 10 4 yrs and ~ 10 5 yrs respectively), and masses 
(<~ 10 M C3 and ~ 10 3 M Q respectively), suggest 
that, if O stars are forming in G20.08N (as it is 
observed), then the smaller scales ought to be re- 
supplied from the larger scales. The same result 
has been found in recent numerical simulations of 
massive star formation in clusters. 



3. The brightest HC H II region (A) has a rising SED 
from cm to mm wavelengths and broad hydrogen 
recombination lines. Both characteristics suggest 
density gradients and supersonic flows inside the 
H II region. A tentative velocity gradient is de- 
tected in the recombination line emission of this 
source, suggesting rotation and outflow in the ion- 
ized gas at the innermost scales. H II region A can 
be interpreted as the inner part of the surrounding 
molecular accretion flow, with the observed ioniza- 
tion being produced by photoevaporation. 



We thank the anonymous referee for a detailed report. 
We also thank Keping Qiu for his comments on the tem- 
perature determination. R. G.-M. acknowledges support 
from an SMA predoctoral fellowship. 



REFERENCES 



Abel, T., Bryan, G.L., & Norman, M.L., 2002, Science, 295, 93 
AfHerbach, A., Churchwell, E., Acord, J. M., Hofner, P., Kurtz, 

S., & Depree, C. G. 1996, ApJS, 106, 423 
Anderson, L. D., & Bania, T. M. 2009, ApJ, 690, 706 
Arce, H. G., Santiago-Garci'a, J., j0rgensen, J. K., Tafalla, M., & 

Bachiller, R. 2008, ApJ, 681, L21 
Avalos, M., Lizano, S., Franco-Hernandez, R., Rodriguez, L. F., 

& Moran, J. M. 2009, ApJ, 690, 1084 
Avalos, M., Lizano, S., Rodriguez, L. F., Franco-Hernandez, R., 

& Moran, J. M. 2006, ApJ, 641, 406 
Beltran, M. T., Cesaroni, R., Neri, R., Codella, C, Furuya, R. S., 

Testi, L., Olmi, L. 2004 ApJ, 601, L187 
Beltran, M. T., Cesaroni, R., Neri, R., Codella, C, Furuya, R. S., 

Testi, L., & Olmi, L. 2005, A&A, 435, 901 
Beltran, M. T., Cesaroni, R., Codella, C, Testi, L., Furuya, R. S., 

& Olmi, L. 2006, Nature, 443, 427 
Bcuthcr, H., Schilkc, P., Gueth, F., McCaughrean, M., Andersen, 

M., Sridharan, T. K., & Menten, K. M. 2002a, A&A, 387, 931 
Beuther, H., Schilkc, P., Sridharan, T. K., Menten, K. M., 

Walmsley, C. M., Wyrowski, F. 2002b, A&A, 383, 892 
Beuther, H., Schilke, P., & Stanke, T. 2003, A&A, 408, 601 
Beuther, H., Zhang, Q., Bergin, E. A., Sridharan, T. K., Hunter, 

T. R., & Leurini, S. 2007, A&A, 468, 1045 
Cesaroni, R., Felli, M., Testi, L., Walmsley, C. M., & Olmi, L. 

1997, A&A, 325, 725 
Churchwell, E. 2002, ARA&A, 40, 27 

Danby, G., Flower, D. R., Valiron, P., Schilke, P., & Walmsley, C. 

M. 1988, MNRAS, 235, 229 
De Pree, C. G., Mehringer, D. M., & Goss, W. M. 1997, ApJ, 

482, 307 

De Pree, C. G., Wilncr, D. J., Mercer, A. J., Davis, L. E., Goss, 

W. M., & Kurtz, S. 2004, ApJ, 600, 286 
Downes, D., Wilson, T. L., Bieging, J., & Wink, J. 1980, A&AS, 

40, 379 

Field, G. B., Blackman, E. G., & Kcto, E. R. 2008, MNRAS, 385, 
181 

Fish, V. L., Reid, M. J., Wilncr, D. J., & Churchwell, E. 2003, 
ApJ, 587, 701 

Franco, J., Kurtz, S., Hofner, P., Testi, L., Garcia-Segura, G., & 

Martos, M. 2000, ApJ, 542, L143 
Galvan-Madrid, R., Rodriguez, L. F., Ho, P. T. P., & Keto, E. 

2008, ApJ, 674, L33 

Garay, G., & Lizano S. 1999, PASP, 111, 1049 

Garay, G., Reid, M. J., & Moran, J. M. 1985, ApJ, 289, 681 

Girart, J. M., Beltran, M. T., Zhang, Q., Rao, R., & Estalella, R. 

2009, Science, 324, 1408 

Gordon, M. A., & Sorochenko, P. N. 2002, Radio Recombination 

Lines: Their Physics and Astronomical Applications 

(Dordrecht: Kluwer) 
Ho, P. T. P., & Haschick, A. D. 1986, ApJ, 304, 501 
Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Ho, P. T. P., Vogel, S. N., Wright, M. C. H., & Haschick, A. D. 

1983, ApJ, 265, 295 
Ho, P. T. P., & Townes C. H. 1983, ARA&A, 21, 239 
Hoarc, M. G., Kurtz, S. E., Lizano, S., Keto, E., & Hofner, P. 

2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, 

& K. Keil (Tucson: Univ. Arizona Press), 181 
Hofner, P., & Churchwell, E. 1996, A&AS, 120, 283 



Hollenbach, D., Johnstone, D., Lizano, S., & Shu, F. 1994, ApJ, 
428, 654 

Jaffe, D. T., & Martin-Pintado, J. 1999, ApJ, 520, 162 
Johnstone, D., Hollenbach, D., & Bally, J., 1998, ApJ, 499, 758 
Kcto, E. 1990, ApJ, 355, 190 
Keto, E. 2002, ApJ, 580, 980 
Kcto, E. 2003, ApJ, 599, 1196 
Kcto, E. 2007, ApJ, 666, 976 

Keto, E., Ho, P. T. P., & Haschick, A. D. 1987a, ApJ, 318, 712 
Kcto, E., Ho, P. T. P., & Haschick, A. D. 1988, ApJ, 324, 920 
Kcto, E., Ho, P. T. P., & Reid, M. J. 1987b, ApJ, 323, L117 
Keto, E., & Klaassen, P. 2008, ApJ, 678, L109 
Keto, E., Welch, W., Reid, M., & Ho, P.T.P. 1995, ApJ, 444, 765 
Keto, E., & Wood, K. 2006, ApJ, 637, 850 
Keto, E., Zhang, Q., & Kurtz, S. 2008, ApJ, 672, 423 
Klaassen, P. D., & Wilson, C. D. 2007, ApJ, 663, 1092 
Klaassen, P. D., & Wilson, C. D. 2008, ApJ, 684, 1273 
Kraemer, K. E., & Jackson, J. M. 1995, ApJ, 439, L9 
Kumar M.S.N., & Grave, J.M.C., 2007, A&A, 472, 155 
Lai, S.-P., Crutcher, R.M., Girart, J.M., & Rao, R., 2001, ApJ, 
561, 864 

Leurini, S., Beuther, H., Schilke, P., Wyrowski, F., Zhang, Q., & 

Menten, K. M. 2007, A&A, 475, 925 
Lizano, S., Canto, J., Garay, G., & Hollenbach, D., 1996, ApJ, 

465, 216 

Lorcn, R. B., & Mundy L. G. 1984, ApJ, 286, 232 

Lugo, J., Lizano, S., & Garay, G. 2004, ApJ, 614, 807 

Mangum, J. G., & Wootten, A. 1994, ApJ, 428, L33 

Mangum, J. G., Wootten, A., & Mundy, L. G. 1992, ApJ, 388, 467 

Martmcz-Garci'a, E. E., Gonzalez- Lopezlira, R. A., & Bruzual-A, 

G., 2009, ApJ, 694, 512 
Molinari, S., Pezzuto, S., Cesaroni, R., Brand, J., Faustini, F., & 

Testi, L. 2008, A&A, 481, 345 
Newman, W. I., & Wasserman, I. 1990, ApJ, 354, 411 
Nummelin, A., Bergman, P., Hjalmarson, A., Friberg, P., Irvine, 

W. M., Millar, T. J., Ohishi, M., & Saito, S. 2000, ApJS, 128, 

213 

Olmi, L., Cesaroni, R., Hofner, P., Kurtz, S., Churchwell, E., & 

Walmsley, C. M. 2003, A&A, 407, 2250 
Patel, N. et al. 2005, Nature, 437, 109 

Plume, R., Jaffe, D. T., & Evans, N. J. 1992, ApJS, 78, 505 
Qiu, K., & Zhang, Q. 2009, ApJ, 702, L66 

Rcmijan, A., Sutton, E. C, Snyder, L. E., Friedel, D. N., Liu, 

S.-Y., & Pei, C.-C. 2004, ApJ, 606, 917 
Rodriguez, L. F., Zapata, L., & Ho. P. T. P. 2009, ApJ, 692, 162 
Roberts, WW., 1969, ApJ, 158, 123 

Scwilo, M., Churchwell, E., Kurtz, S., Goss, W. M., & Hofner, P. 

2004, ApJ, 605, 285 
Sewilo, M., Churchwell, E., Kurtz, S., Goss, W. M., & Hofner, P. 

2008, ApJ, 681, 350 
Shu, F.H., Milione, V., Gebel, W., Yuan, C, Goldsmith, D.W., & 

Roberts, WW., 1972, ApJ, 173, 557 
Shu, F. H., Adams, F. C, & Lizano, S. 1987, ARA&A, 25, 23 
Sollins, P. K. 2005, PhD thesis, Harvard University, ISBN: 

3-936586-35-7 

Sollins, P. K., Zhang, Q., Keto, E., & Ho, P. T. P. 2005a, ApJ, 
624, L49 



10 



GALVAN-MADRID ET AL. 



Soiling, P. K., Zhang, Q., Keto, E., & Ho, P. T. P. 2005b, ApJ, 
631 399 

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 
319 730 

Sridharan, T. K., Beuther, H., Schilke, P., Menten, K. M., & 

Wyrowski, F. 2002, ApJ, 566, 931 
Tang, Y.-W., Ho, P.T.P., Koch, P.M., Girart, J.M., Lai, S.-P., & 

Rao, R., 2009, ApJ, 700, 251 
Turner, B. E. 1979, A&AS, 37, 1 

Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 
Vazquez-Semadeni, E., Gomez, G.C., Jappsen, A.K., 

Ballesteros-Paredes, J., & Klessen, R.S., 2009, arXiv:0904.4515 
Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 

1998, MNRAS, 301, 640 
Walsh, A. J., Longmore, S. N., Thorwirth, S., Urquhart, J. S., & 

Purcell, C. R. 2007, MNRAS, 382, L35 
Wang , P., Li, Z.-Y ., Abel, T., & Nakamura, F. 2009, 

arXiv:0908.4129 



Wilner, D. J., Wright, M. C. H., & Plambeck, R. L. 1994, ApJ, 
422 642 

Wood', D. O. S., & Churchwell, E. 1989, ApJS, 69, 831 

Wu, J., & Evans, N. J., II 2003, ApJ, 592, L79 

Young, L. M., Keto, E., & Ho, P. T. P. 1998, ApJ, 507, 270 

Yorke, H. & Welz, A., 1996, A&A, 315, 555 

Zhang, Q., & Ho, P. T. P. 1995, ApJ, 450, L63 

Zhang, Q., & Ho, P. T. P. 1997, ApJ, 488, 241 

Zhang, Q., Ho, P. T. P., & Ohashi, N. 1998, ApJ, 494, 636 

Zhang, Q., Hunter, T. R., Beuther, H., Sridharan, T. K., Liu, 

S.-Y., Su, Y.-N., Chen, H.-R., & Chen, Y. 2007, ApJ, 658, 1152 
Zhang, Q., Hunter, T. R., Sridharan, T. K., & Ho, P. T. P. 2002, 

ApJ, 566, 982 



HIERARCHICAL ACCRETION IN G20.08-0.14 N 



TABLE 1 

Observational Parameters 



Epoch 


Array 


Phase Center 3, 


Bandpass 


Phase 




Flux 








a(J2000) 


d(J2000) 


Calibrator 


Calibrator 


Calibrator 


2003 Apr 28 


VLA-D 


18 28 10.384 


-11 28 48.65 


3C454.3 


1833 - 


210 


0137 + 


331 


2003 May 13 


VLA-D 


18 28 10.384 


-11 28 48.65 


3C454.3 


1851 + 


005 


1331 + 


305 


2003 Oct 09 


VLA-BnA 


18 28 10.384 


-11 28 48.65 


3C454.3 


1851 + 


005 


1331 + 


305 


2003 Oct 10 


VLA-BnA 


18 28 10.384 


-11 28 48.65 


3C454.3 


1851 + 


005 


1331 + 


305 


2006 Jun 25 


SMA-VEX 


18 28 10.38 


-11 28 48.60 


3C273 


1830 + 


063 


1830 + 


063 


2006 Jun 25 


SMA-VEX 


18 28 10.76 


-11 29 27.60 


3C273 


1830 + 


063 


1830 + 


063 


2006 Jul 06 


SMA-VEX 


18 28 10.38 


-11 28 48.60 


3C454.3 


1751 + 


096 


1751 + 


096 


2006 Jul 06 


SMA-VEX 


18 28 10.76 


-11 29 27.60 


3C454.3 


1751 + 


096 


1751 + 


096 


2007 Oct 26 


VLA-B 


18 28 10.400 


-11 28 49.00 


3C454.3 


1743 - 


038 


1331 + 


305 


2007 Oct 27 


VLA-B 


18 28 10.400 


-11 28 49.00 


3C454.3 


1743 - 


038 


1331 + 


305 



a Units of right ascension are hours, minutes, and seconds. Units of declination arc degrees, arcminutcs, 
and arcseconds. 



TABLE 2 

LlNES a 



Species 


Transition 


^0 


Array 


HPBW 






(GHz) 




(arcsec X arcsec; deg) 


H 


66a 


22.364178 


VLA-B 


0.47 X 0.34; 3 


NH 3 


(2,2) 


23.722633 


VLA-D 


4.68 x 2.93; 236 


NH 3 


(3,3) 


23.870129 


VLA-D 


4.50 X 3.32; 8 


NH 3 


(2,2) 


23.722633 


VLA-BnA 


0.37 x 0.28; 343 


NH 3 


(3,3) 


23.870129 


VLA-BnA 


0.71 X 0.34; 


13 CO 


2-1 


220.398681 


SMA-VEX 


0.55 X 0.41; 32 


CH 3 CN 


12 (4)- 11 (4) 


220.679297 


SMA-VEX 


0.55 X 0.41; 32 


CH 3 CN 


12(3)- 11(3) 


220.709024 


SMA-VEX 


0.55 X 0.41; 32 


CH 3 CN 


12(2)-11(2) 


220.730266 


SMA-VEX 


0.55 X 0.41; 32 


S0 2 


11(1,11)-10(0,10) 


221.965200 


SMA-VEX 


0.54 X 0.41; 32 


CO 


2-1 


230.538000 


SMA-VEX 


0.53 x 0.39; 37 


OCS 


19-18 


231.060991 


SMA-VEX 


0.53 x 0.39; 37 


H 


30a 


231.9009 


SMA-VEX 


0.53 x 0.39; 37 



CH 3 CN 
HNCO 
CH 3 CN 
CH 3 CN 
CH 3 CN 
CH 3 CN 
CH 2 CHCN 
13 CS 



12(7)-11(7) 
10(1,9)-9(1,8) 
12(6)- 11(6) 
12(5)-11(5) 
12(l)-H(l) 
12(0)-11(0) 
24(0,24)-23(0,23) 
5-4 



220.539340 
220.584762 
220.594438 
220.641096 
220.743015 
220.747265 
221.76598 
231.220768 



SMA-VEX 
SMA-VEX 
SMA-VEX 
SMA-VEX 
SMA-VEX 
SMA-VEX 
SMA-VEX 
SMA-VEX 



0.55 X 0.41; 32 
0.54 X 0.41; 32 
0.55 X 0.41; 32 
0.55 X 0.41; 32 
0.55 X 0.41; 32 
0.55 X 0.41; 32 
0.55 X 0.41; 32 
0.52 x 0.39; 37 



a Detected lines at S/N > 6. The top part of the table lists the lines with S/N > 10 
(except the H66q line) and clearly isolated in frequency. The bottom part of the table lists 
the lines detected at S/N < 10 or blended. Some spectral features (see Fig. that were 
not properly identified due to low S/N (~ 5) and blending are not listed. 
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TABLE 3 

Model for the H II Region G20.08N A 



Paramctcr a 


Value 


HII Radius (AU) 


2530 


Electron density 1 " (10 5 cm -3 ) 


1.4 


Exponcnt c 


1.3 


Gas Mass d (M ) 


35-95 


Spectral type 


07.5 


Stellar Mass (M ) 


34 


HII Mass (M Q ) 


0.05 



a The independent parameters are the first 
four rows. 

Electron density at the HII radius. 
c Exponent 7 of the power-law density gradi- 
ent in the ionized gas, where n oc r~ " y . 
^ Mass of molecular gas obtained from the 
dust emission. Assuming dust temperature 
Trf — 230 K and a gas-to-dust ratio of 100. 
The range is caused by the dust emissivity co- 
efficient (3 used (see IKeto et aLl bOOgl). from 
/3 = 1 to 13 = 1.5. 



TABLE 4 

Emission Line Parameters Toward H II Region A a 



Species 


Transition 


Vlsr 


FWHM 


-^B ,peak 






( km s" 1 ) 


( km s' 1 ) 


(K) 


H 


66« 


39.3 ± 1.1 


34.1 ± 2.5 


190 ± 13 


H 


30a 


40.3 ±0.4 


32.5 ±0.9 


89 ±2 


CH3CN 


12(4)-11(4) 


42.7 ±0.3 


5.8 ±0.7 


32 ±4 


CH3CN 


12(3)-11(3) 


41.4 ±0.2 


5.3 ±0.5 


41 ±4 


CH3CN 


12(2)-11(2) 


41.4 ±0.2 


4.9 ±0.5 


32 ±3 



a From Gaussian fits, lo" statistical errors are quoted. 
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231.0 231.5 
Rest Frequency [GHz] 




221.0 221.5 
Rest Frequency [GHz] 



Fig. 1. — Wide-band, continuum-free spectra from the SMA-VEX data at the position of the 1.3-mm continuum peak (see Table 2 for 
details). The channel spacing in this plot is 3 km s . The question mark (?) in the upper sideband (top frame) might be a superposition 
of lines of CH3OCH3 and CH3CH2CN. The ?s at the center of the lower sideband (bottom frame) could be from vibrationally excited 
CH3CN. The ? in the upper sideband close to the SO2 might be from CH3CCH. 



o 




18 h 28 m 10 s 55 10 s 50 10 s 45 10 s 40 10 s 35 10 s 30 10 s 25 
RA (J2000) 



Fig. 2. — VLA-BnA 1.3-cm continuum (red contours) overlaid on the SMA-VEX 1.3-mm continuum (color scale) toward the G20.08N 
complex. From west to east: H II region A is the compact, strongest peak at both wavelengths. H II region B is the less bright H II 
region at < 1" to the SE of A. H II region C is the more extended emission to the SE of B. The color scale goes linearly from 8 to 430 
mjy beam -1 (the rms noise in the mm image is 2 mjy beam -1 ). Contours are placed at —5, 5, 10, 20, 40, 80, 150 X 1 mjy beam -1 , the 
noise of the cm image. The SMA-VEX beam (empty ellipse) encircles the VLA-BnA beam (filled ellipse) at the bottom-left of the image. 
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9.0 9.5 10.0 10.5 11.0 
log Frequency (Hz) 



11.5 12.0 



Fig. 3. — Radio-to-millimctcr SED of G20.08N A. The 6-cm and 2-cm measurements were obtained from Wood & Churchwcll (1989). 
The 1.3-cm and 1.3-mm points were obtained from Gaussian fits to our VLA-BnA and SMA-VEX data, respectively. Crosses are the data 
points. The error bars correspond to the 10% and 20% uncertainty expected in the VLA and SMA flux measurements, respectively. The 
red dashed line shows the flux of an H II region with a density gradient. The black dotted line is the flux from the warm dust component. 
The solid blue line is the sum of the two components. The relative contributions of free-free (70 %) and dust (30 %) to the 1.3-mm flux 
were estimated from the observed H30o line-to-continuum ratio, and are marked with diamonds. 
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Fig. 4. — Channel maps of the VLA-D NH3 (3,3) observations. Emission is in solid contours and absorption in dashed contours. The 
star covers the H II regions shown in Fig. [T] Contour levels are at —35, —25, —15, —10, —7, —5, 5, 7, 10, 15, 25, 35 X 2 mjy beam -1 . A 
clear velocity gradient in emission is seen from one side of the absorption to the other. The LSR systemic velocity of the molecular gas is 
V S y a = 42.0 km s _1 . The original maps at 0.3 km s 1 spectral resolution were smoothed to 0.6 km s 1 for clarity. 
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Fig. 5. — Position-velocity diagrams of NH3 (2,2) and (3,3) from the VLA-D data. The cuts were done at PA=45° 
(bottom row) and PA=135° (top row). Dashed contours are absorption, solid contours are emission. Contouring is at 
-168, -144, -120, -96, -72, -48, -36, -24, -12, -4, 4, 8, 12, 16, 24, 32, 40, 48, 56, 64 X 1 mjy beam -1 . Only one inner satellite line is cov- 
ered in the velocity range. In the SW-NE cuts (across the major axis of the cloud) the difference in the velocity of the emission with respect 
to V 3 ys = 42 km s — 1 increases closer to the position center. This can be interpreted as spin up with decreasing distance from the center. 
However, the same trend is present in the NW-SE cuts (along the minor axis of the cloud), although only in the stronger, blueshifted side 
of the emission. This suggests that besides rotation, radial motions in the frame of the central stars are also present. 





E _5 
« 

-15 
-20 

15 20 25 30 35 40 45 50 
V LSR Ikm/s] 





g -5 

I -10 

■d 
« 

-15 
-20 

15 20 25 30 35 40 45 50 
Vlsr [km/s] 

Fig. 6. — NH3 spectra from the VLA-D observations toward the center of the absorption in G20.08N. The vertical line marks the 
systemic velocity (V sy s = 42 km s -1 ). The small peak near the middle of the (3,3) absorption is due to maser emission (confirmed in 
the high-resolution, VLA-BnA data, see £|3.6| I. The absorption peaks in the main (2,2) and (3,3) lines are redshifted with respect to V sy s, 
indicating the presence of inflow in the kinematics of the parsec-scale molecular cloud. The other weak absorption component seen in the 
(2,2) spectrum is one of the inner satellites. The absorption component at rj 22 km s -1 in the (3,3) spectrum is also an inner satellite. 
The weaker absorption at 30 — 37 km s" 1 in the (3,3) spectrum may arise from an outflow. The blueshifted absorption in the main lines 
closest to Vsys might arise from blending with the rotation seen in emission. 
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Fig. 7.— Channel maps of OCS J = 19 - 18 from the SMA-VEX observations. Contours are -4,4,6,8,10,12,14 X 30 mjy beam" 1 
(negative in dashed, and positive in solid). The peak intensity is 423 mjy beam" 1 . The numbers in the upper left corner indicate the 
central LSR velocity of the channel. The two stars mark the positions of H II regions A (west) and B (east). 
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Fig. 7 (continued). — Channel maps of SO2 J{K a ,K b ) = 11(1,11) - 10(0,10) from the SMA-VEX observations. Contours are 
—4, 4, 6, 8, 10, 12, 14, 16, 18, 20 X 30 mjy beam -1 (negative in dashed, and positive in solid). The peak intensity is 612 mjy beam -1 . 
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Fig. 7 (continued).— Channel maps of CH 3 CN J(K) = 12(3)-11(3) from the SMA-VEX observations. Contours are -4, 4, 6, 8, 10, 12, 14 
X 30 mjy beam -1 (negative in dashed, and positive in solid). The peak intensity is 440 mjy beam -1 . 
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Fig. 7 (continued).— Channel maps of CH 3 CN J(K) = 12(4)-11(4) from the SMA-VEX observations. Contours are -4, 4, 6, 8, 10, 12, 13 
X 30 mjy beam -1 (negative in dashed, and positive in solid). The peak intensity is 401 mjy beam -1 . 
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Fig. 8. — Velocity-integrated emission (moment 0, contours) and intensity-weighted mean velocity (moment 1, color scale) maps from 
hot-core molecules toward G20.08N. The Top Left panel shows the OCS J = 19 — 18 (contours are —5, 5, 9, 13, 17 X 0.15 Jy beam -1 km s -1 ). 
The Top Right panel corresponds to the S0 2 J(K a , K b ) = 11(1, 11) - 10(0, 10) (contours are -5, 5, 8, 12, 15 X 0.25 Jy beam -1 km s -1 ). 
The Bottom Left panel shows the CH 3 CN J(K) = 12(3) -11(3) (contours are -5,5,8, 12, 15 X 0.15 Jy beam -1 km s -1 ). The Bottom Right 
panel plots the CH3CN J(K) = 12(4) — 11(4) (contours are —5,5,9,14,19 X 0.1 Jy beam -1 km s -1 ). Negative contours are in dashed 
style, and positive contours in solid. Only the continuum H II region A is associated with warm gas in our SMA-VEX data. The extent 
of the line emission is similar in all the tracers, although for the SO2 it is more extended toward the northwest and more absorbed toward 
the continuum peak. The velocity gradient seen in the channel maps for each molecule is also seen here. The color scale is the same in all 
the frames. 
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Fig. 9. — Velocity-integrated emission (moment 0, contours) and velocity dispersion with respect to moment-1 velocity (moment 2, 
gray scale) maps from hot-core molecules toward G20.08N. Contours are the same as in Fig. [8] The plotted velocity dispersion is 
cr = FWHM/2(2 ln2) 1 / 2 . a increases toward the center of H II region A, possibly caused by unresolved motions toward the center. 
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Fig. 10. — Position-velocity diagrams of the SMA-VEX observations. The left column is for a cut at PA= 45° (SW-NE). The right 
column is at PA= 135° (NW-SE). Cuts are centered at the position of the continuum peak of H II region A. Contours are the same as 
in Fig. [7] The velocity gradient interpreted as rotation from SW to NE is clearly seen in OCS and CH3CN. The SO2 features are more 
complex. 
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Fig. 11. — Spectra (points) and Gaussian fits (lines) to the CH3CN J = 12 — 11, if = 2, 3, and 4 lines at the position of H II region A. 
Error bars denote the la noise in the 2 km s _1 -wide channels. 
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Fig. 12. — Top and middle panels: NH3 (2,2) and (3,3) spectra toward the absorption peak of H II region A in the VLA-BnA observations. 
The vertical line marks the systemic velocity (V ay s = 41.8 km s^ 1 ) of the molecular gas at scales comparable to the H II region. The 
centers of the absorption lines are redshifted with respect to V sya , indicating inflow of molecular gas toward H II region A at small scales. 
The (3,3) absorption is broader (FWHM 6.5 km s —1 ) than the (2,2) (FWHM 5.6 km s — 1 ), probably caused by larger motions closer 
to the center. Bottom panel: NH3 (2,2) spectrum toward H II region B. The absorption, considerably fainter than for H II region A, is 
also redshifted. 
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Fig. 13. — Rotation diagram for the CH3CN J = 12 — 11. The top panel is a fit of all the K = 0, 7 lines. The bottom panel uses only 
the i<" = 4, 7 lines, which have lower optical depths. Error bars are 3cr values. Lines are the linear fits to the data. 
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Fig. 14. — Velocity-integrated emission (moment 0, contours) and intensity-weighted mean velocity (moment 1, color scale) map of the 
H30a RRL emission toward G20.08N. Contours are at 7, 50 X 0.5 Jy beam" 1 km s _1 . Although the emission is unresolved at half power, 
the moment 1 map hints at the presence of a velocity gradient in the ionized gas similar to that seen in CH3CN and OCS. 
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Fig. 15. — Spectra (points) and Gaussian fits (lines) to the H30a (top panel) and H66« (bottom panel) lines toward G20.08N A. Error 
bars denote the lc noise in the channels. The channel spacing is 2 km s for H30a and 5.2 km s — 1 for H66a. The flux was integrated 
over a 0.5" square box centered on H II region A. 
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Fig. 16. — VLA-BnA spectrum of the NH3 (3,3) maser spot toward G20.08N. The brightness temperature Tb scale assumes an angular 
size equal to half the uniform-weighted beam dimensions. 



